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Summary 



This report gives briefly the theories that form the basis for calcu- 
lation of the design and the behaviour of a windmill. A modification 
of the Prandtl model of tip losses is derived. This modification takes 
the relatively heavy loading of the windmill rotor into account. It is 
argumented that, in contrast with propeller design, a maximum energy 
extraction is reached by enlarging the chords of the blades near the 
tips. Selection, design and construction of several roitors and of a 
test unit are described. Tests of steel plate rotors and so-called 
sail trouser rotors are described while test results are presented 
in the form of C p -A and Cg-X characteristics. Final conclusion is 
that with simple materials high power coefficients are possible. 
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1. THEORY 



In the past few years several authors translated the propel ler- theory, as 
developed before World War II, into theories that can be used for design, 
calculation and prediction of the behaviour of windmills (l, 5, 8). In 
this chapter I have tried to expose the line of reasoning and the essen- 
tial results from the above given references. Special consideration has 
been given to the effect of a finite number of blades, since the existing 
models to describe this effect were not found satisfactory. 

1.1. Momentum theory 



1.1.1. Axial momentum theory (l) 

The simplest description of the extraction of energy from the wind is the 
one dimensional, incompressible, non-viscous flow model using axial mo- 
mentum theory (Rankine-Froude) (see fig. 1.1.). The flow is assumed to be 
entirely axial with no rotational motion. Two expressions for the thrust 
acting upon the extracting device, may be obtained: 



Continuity 
Momentum theorem 
Bernoulli upwind 
Bernoulli downwind 
thus p - p' = iP (V, 



UA = 

T - m (V.- U x ) = pAU (V w - U x ) 

P + iP yj = P + iP U 

p + Jp u\ 2 = p' + iP u 2 

2 2 

. - u, ) 



1-la) 

1-1) 

1-2) 

;i-3) 
:i-4) 



Fig. 1.1.: Axial flow model 
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Combining equations (1-1) and (1-4) with: 
T = AAp 



leads to: U = 



V + U 

OC j 



The power extracted is: 
= IP UA (vj - u/) 

F 

If we denote C n = 



and U = (1-a) V ro then with (1-6) = (l-2a) V 
and with (1-7) C p = 4a (1-a) 2 (see fig. 1.1. a.) 
which has a maximum when a = ~ 

C Prnax = TI ' °' 593 



(1-5) 
(1-6) 



(1-7) 
(1-8) 

(1-9) 




Fig. 1.1. a. : 

The factor (a) is known as the axial interference factor and is a measure 
of the influence of the turbine on the speed of the air. Note that even if 
we don't know with which device we shall extract the energy, we al ready do 
know that the highest possible energy extraction requires slowing down the 
air to one third of its original speed. 

16 

In the used model, that gives the maximum Cp equal to ^y, three principle 
assumptions have been made: 
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1. The flow is entirely axial, 

2. no friction occurs when the air passes the device, 

3. the flow is rotational ly symmetric. 

If the device that extracts power from the air is a rotor, the power of 
the air is transformed into torque Q and angular speed ft. This means that, 
if the angular speed cannot be infinite, extracting power {= P = Q * ft) 
means generating torque. 

Since generating torque implies tangential forces and thus momentum changes 
in the air in tangential direction, assumption 1 is only acceptable for 
very high speed devices. The higher the torque that is generated, the higher 
the tangential momentum in the air downstream will be. Apart from other los- 
ses that occur in reality, we here have the first reason why a Cp of 59.3% 
cannot be realized in a real construction. Deviations in reality from the 
second and third assumption will lead to even more losses. These losses will 
be iealt with in the following paragraphs. 

1.1.2. Effect of wake rotation 

In describing this effect the assumption is made that upstream of the rotor 
the flow is entirely axial and that the flow downstream rotates with an 
angular velocity w (r), but remains irrotational . 

According to (8) this model produces unrealistic velocities near the rota- 
tion axis; these can be avoided by replacing the region of high angular ve- 
locities by a rigid rotating core. In the following discussion we will not 
include this effect because it will not affect the results significantly. 




Fig. 1.2.: Axial flow model with wake rotation. 
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Orie of the results of the axial flow model was that: - 
V + U 

U = " * U * - (1-6) 
I 

According to (8) this is not correct for the model with wake-rotation; a 
graph presented in [8] indicates that holding on to the relation (1-6) 
will not result in too large an error for practical cases (X > 1.5) (see 
Appendix I.). 

Expressions for torque and maximum power may now be obtained by considering 




Fig. 1.3. ; Annul us at radius r. 

the flow through an annul us at radius r with area AA = 2-rrrAr. By changing 

the momentum in the air in tangential direction, tangential forces act 
upon the rotor: 



AF t = m (AV) = pUAAwr 

= 27rpUojr 2 A r (1-10) 
The torque generated in the annul us Ar is: 

AQ = 27rpUwr 3 Ar (1-11) 

And since power = torque * angular speed, the power extracted is: 

AP = 2irpfiUu3r 3 Ar (1-12) 

The torque of the rotor becomes: 
R 

Q = 2irp / Our dr (1-13) 


and the rotor power: 
R 

P = 27rpft / Uor dr (1-14) 
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To find the maximum of the expression (1-14) a relation between U and w 
is needed. This relation can be found by equating the axial force on an an- 
nul us at radius r as given by the axial momentum theory to the axial force 
on the same annul us as may be found from considering the circulation at ra- 
dius r. This is done in fl] and the result is: 

U (V^ - U) = icor (£2 + |o))r (1-15) 

On page 13 this result will be found in another way in dimensionless form. 
Defining a and a' by: 

U = (1-a) (1-16) 

ur = 2a'ftr (1-17) 
and 

A r = fir/V M (1-18) 

A = fiR/V ro (1-19) 

C p = (1-20) 

00 

We may obtain obtain from (1-14) and (1-15): 

C D = / (1-a) a'A 3 dA (1-21) 
A O r r 

and (l-a)a = a' (1 + a')X * (1-22) 

Now the power integral of (1-21) can be made stationary subject to the con- 
straint of (1-22); see (l, 8] . The result is a relation between a and a'. 

a' = (l-3a)/(4a-l) (1-23) 

This result will be used in paragraph 1.3 for design purposes. 

Hence ^ < a < ~ . By varying a between these extremes in equations (1-13), 
(1-22) and (1-21), C p for various values of A may be obtained. The results 
are shown in table 1.1. 
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a 


a' 


X 

r 




X 


C P 


Cpd-24) 


0.25 


GO 


0.000 







0.0 


0.0 


0.26 


5.5 


0.073 




0.5 


0.288 


0.249 


0.27 


2.375 


0.157 




1.0 


0.416 


0.416 


0.28 


1.333 


0.255 




1.5 


0.481 


0.481 


0.29 


0.812 


0.374 




2.0 


0.513 


0.513 


0.30 


0.500 


0.529 




2.5 


0.533 


0.532 


0.31 


0.292 


0.753 




5.0 


0.570 


0.567 


0.32 


0.143 


1.154 




7.5 


0.582 


0.577 


0.33 


0.031 


2.619 




10.0 


0.585 


0.582 


1/3 


0.000 






CO 


0.592 


0.592 



Table 1.1. : a) relation between a, a' and X^ for ideal flow and b) Cp = 
f{X) for ideal flow. 

For the present work it was useful to derive an approximation of Cp, as 
given by th* formula: 



16 -0.3538X 



-1.2946 



'P approx 27 



(1-24) 



As can be seen from the table 1.1 this approximation fits rather well for 
X <t 1 (error < 1 percent). 

Fig. 1.4. shows C p = f(X) in graphical form. Keeping in mind the two assump- 
tions in the model (rotationally symmetric and frictionless flow) this Cp 
may be interpreted as the characteristic of an ideal windmill. 




1 2 3 4 5 
Fi 9- M - : C P-n,ax as f < X >- 
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As will be seen in the next chapters, a windmill would need, apart from a 
frictionless medium and an infinite number of blades, an unlimited possibi- 
lity of changing the blade angles and chords if a characteristic, as in fig. 
1.4. , were to be realized. Since this is not possible, when thinking of the 
horizontal-axis windmills we have in mind, it is, in contrast to the inter- 
pretation often found in other literature, better to interpret fig. 1.4. 
as the collection of maximum C p -values of the characteristics of different 
windmills, each operating at its specific optimum X. 

Blade element theory 

Since the device that extracts the energy in the momentum theory is not spe- 
cified, the momentum theory cannot give information on: 

- blade chords and angles; 

- influence of friction; 

- influence of number of blades. 

The method that, when used in conjunction with the momentum theory, gives in- 
formation about these aspects is called BLADE ELEMENT THEORY. The basic as- 
sumption is that the flow through the rotor occurs in non-interacting circu- 
lar stream tubes. 

The flow at an element Ar (see fig. 1.5.) may then be regarded as two dimen- 
sional. This opens the possibility of using the information that exists on 
the properties of airfoil sections, since that information is in fact strict- 
ly for two-dimensional flow. 



Fig. 1.5. : Rotor blade elements. 

In this paragraph the flow through the circular stream tube is thought to be 
rotational ly symmetric. Since the number of blades is always finite this as- 
sumption is not entirely valid and the formulas found here will be corrected 



/ 




R 



-13- 



in the following paragraphs. 

1.2.1. Effect of drag 

An element at radius r in fig, 1.5. is shown in fig. 1.6. In this figure 
the rotational velocity is assumed to be half the value of the final rota- 
tional velocity in the wake. 



Defining: 



C D = D/|pw 2 A 


(1-25) 


C L = L/£pW 2 A 


(1-26) 


C x = F x /|pW 2 A 


(1-27) 


C y = F y /£pW 2 A 


(1-28) 


wi th : 





A = cAr and c = blade chord, 




Fig. 1.6. : Velocities and forces on blade element at radius r. 

the following trigonometric relations may be obtained from fig. 1.6.: 

a = <J> - B (1-29) 

tan* = ) ' a , ^ (1-30) 

1 + a Qr 

Cy = C^cos<f> + Opsins (1-31) 

C x = C^sin* - C D cos0 (1-32) 
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Expressions may now be derived for the thrust T and torque Q of element Ar at 
radius r: 

dT - B c JpvfCydr (1-33) 
dQ = B c £pW 2 C x r dr (1-34) 

where B is the number of blades. From the momentum theory (1-1) and (1-11) 
we already know: 

dT = p(2TTrdr) U (V^ - UJ (1-1) 
dQ = p(2iTrdr) Owr* (1-11) 

Equation (1-1) and (1-11) were derived from the momentum theory. Since this 
theory essentially does not include friction, we may only equate the thrust 
and torque of equation (1-1) and (1-11) to equation (1-33) and (1-34) if we 
assume = 0. 

Using definitions (1-16) and (1-17) the following equations may be derived: 

aC, cosd) 

-W- (1-35) 



1 " a a *;- z 

4 Sin 



a' aC L 



1 + a' 4 cos<f) 
Be 



(1-36) 



with 



c 



27Tr 



Note that taking together (1-35) and (1-36) leads with (1-30) to equation 
(1-15) in dimensionless form of (1-22). From (8) the following elegant me- 
thod for calculating the flow conditions at an element Ar at radius r is 
given. 

Given r, c, CJa) , C D (a), S, V w , fi, B: 

A. Estimate a and a' 

B. Calculate <j> (1-30) 

C. Calculate a (1-29) 

D. Find Cj_ 

E. Calculate a (1-35) 

F. Calculate a' (1-36) 

G. Go back to step B and repeat. 



Once this iteration converges, the sectional flow properties are known and 
the local contributions to torque (1-34) and thrust (1-33) may be integrated 
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to determine the overall torque and thrust of the rotor. Note that now dQ 
in equation (1-34) must be interpreted as net torque and that for we 
need equation (1-32) as a whole. 

The effect of the profile drag on the C p - X curve of the ideal windmill 
(see page 10) can be found through the effect on the torque (l). 



Without drag dQ = BcJpW C^sin<f>rdr 



with drag dQ - Bc|pW (C^si n<f> - C^cos*) rdr 

C ■ 



(dQ) 



with drag 



= Bc£pW C L sin<f> 

1 - jr~ COt(|) 



1 - 



-J^ COt(J) 



rdr 



without drag 



equation (1-21) becomes: 

X C 
C p = % / (1 - a)a'X 3 (1 - D cot*) dX = 
X r \ r 



= c 



P ideal ^"D 

X 



(C n = o) - c 



where C w = ~ / (1 - a)a'X " J± cot* 

w X r u r 



W 



dX. 



(1-37) 



together with (1-30) C w becomes: 

C w =5_ /(l-a)a'X r 3 ^ |±f X p dX r 

Assuming that the angle of incidence at the blade section and the form of the 
section itself do not vary along the blade. 



W ,2 



8 C D X 



^ / a'(l + a')X d X 
\ r r 



with equation (1-22) 



w X L L r 1 



With table (1.1.) it may be concluded that for X^ > 1 under the optimum con- 
ditions (1 - a)a is almost constant (error < 0.25%) and equal to 0.22. 



now becomes: 

c ^ 8 C D ) 2 A a . , 
X L L y r r 



* 16 X 

* 77 A 



(1-38) 
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In fig. 1.7. Cp - X curves are shown for different Cp/C^-values. The curves 
should be interpreted in the way that is indicated on page 11. 

.2. Effect of a finite number of blades 

The theories presented in the previous paragraphs, are based on the assump- 
tion that the flow through the rotor is rotationally symmetric. With a finite 
number of blades this assumption obviously does not hold. Another assumption 
is that the flow at any element is two-dimensional. 

Radial acceleration and wake-induced flow at the tip may alter the assumed 
flow pattern. According to (8), the radial acceleration may be neglected 
for most wind power machines; however, the wake effects may not be neglec- 
ted. Since the effects, due to the finite number of blades, result in per- 
formance losses concentrated near the tips of the blades, we speak of tip 
losses and we describe these effects with tip loss models. 
The tip losses have been analysed by Prandtl (3). He developed a model des- 
cribing the reduction of circulation due to the wake interaction at the tips. 
This model was derived for lightly loaded propellers with neglegible wake 
contraction. "Lightly loaded" should be interpreted as "giving a small change 
in the momentum of the flow that passes through the rotor". 
The windmill is destined to produce as much power as possible. We know, from 
the momentum theory, that we therefore must reduce the windspeed to about one 
third of its original value. This means that the assumption of light loading 
may lead to erroneous results in calculating windmill performance. 

Although this means that one of the basic assumptions of Prandtl's analysis 
may not be valid for the windmill, we will use here his results because it 
will be shown that some changes can be made that take the relatively heavy 
loading into account. 

The tip losses are expressed by a circulation-reduction factor defined by: 
F -{£ (1-39) 

00 

where r is the actual circulation of one blade at radius r, B is the number 
of blades and r is the circulation of a rotor with an infinite number of 
blades as calculated by the previously mentioned theory. The factor F is a 
function of the number of blades, tip speed ratio and radial position. It has 
been derived for a frictionless rotor with optimum distribution of circula- 
tion along the blades. 
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A more refined model was given by Goldstein, but we will use the Frandtl 
model for the following reasons: 

a) The difference in results of the two models is small (8); 

b) The Goldstein model is very difficult to use while the Prandtl model is 
relatively simple; 

c) The above given considerations on the assumptions with regard to the 
load apply to both models. 

The Prandtl model is described in (3, 5). The result is a formula for factor 
F: 

F = | arccos e" f (1-40) 
where 

f = B (1 _ £ )N /TT? (1-41) 

The quantity\/l + X* came i nto Prandtl 's model when he needed the distance 
between the vortex sheets in the wake of the flow. See fig. 1.8., where for 
a lightly loaded one bladed propeller, the vortex sheet is shown. 




Fig. 1.8. : Vortex sheet of lightly loaded one bladed propeller. 

For more blades the situation is shown in fig. 1.9a. where the distance s 
between the sheets is: 

s = ™ sin* 

where <f> is the angle of the vortex sheets (fig. 1.9a.) which may be inter- 
preted as the angle <j> in fig. 1.6. 
In the non-contracting wake of the Prandtl model: 

v 1 

oo 1 
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Fig. 1.9a. : Vortex sheets of lightly loaded multi bladed propeller. 

Since the Prandtl model considers the wake of the flow we think it will be 
the best to take for <j> the angle that occurs in the wake of the not lightly 
loaded rotor (see fig. 1.9b.). Analogous to equation (1-30) this angle ^ 
may be calculated from: 

tan$l = f^f, X L (1-30.) 

According to (3, 8) the physical meaning of the tipspeed, correction is virtu- 
ally that the maximum change of axial velocity, (V^ - or 2aV oo , in the 
slipstream only occurs on the vortex sheets and the average velocity change 
is only a fraction F of this velocity. Thus the axial velocity change 2aV oo 
becomes 2aFV and in a similar manner, the angular velocity change is writ- 

CO 

ten 2a'Ffi. Note that for calculating velocity changes on the blade element 
itself the velocity change remains 2aV oo and 2a'ft but for calculating momen- 
tums we should use 2aFV and 2a'Fft. 

00 

Therefore (1-33) and (1-34) remain the same but in equating them to (1-1) and 
(1-11) we should take in (1-1) and (1-11) for U 1 and 03: 

U 1 = (1 - 2aF) (1-44) 

03 = 2a'FS7 (1-45) 

In (8) it is suggested that the correction of the velocity distribution should 
not only be applied to the velocity u\ of the wake, but also to the velocity 
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U in the plane of the rotor, for calculating the average mass flow through 
an annul us. 
U = (1-aF) V ro 




fig. 1.9b. Vortex sheets of heavy loaded multy bladed propeller. 

Equating again (1-1) and (1-11) to (1-33) and (1-34) changes the relations 
(1-35) and (1-36): they now become 



aF (1 - aF) = 



aC L cos<t> (1-a)' 
4s in^ cf> 



(1 - 47) 



and a ' F d - * F > 
(1 - a)(l + a') 



aC, 



4 cos* 



(1 - 48) 



The method for calculating flow conditions and torque contributions, as 

described on page 13, only changes in that for steps E and F now equations 
(1 - 47) and (1 - 48) should be used. 

Factor F is to be calculated from (1 - 40) and (1 - 43). 
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Note, that in combining (1 - 47) and (1 - 48) again equation (1 - 22) may 
be derived. This means that the Prandtl factor F does not change the 
relation between interference factors a and a'. In trying to find a 
formula for the effect of the finite number of blades on the maximum 
power of the propeller as a whole, Prandtl derived an effective radius for 
lightly loaded optimum propellers. If the equations as derived from 
momentum theory are used for calculating form and angles of the blades 
of the windmill, the above mentioned correction may be used for calculating 
the maximum power. The effective radius may be calculated with 

Re = x . 1,386 sin ^ (1 _ 4g) 



Since the total power increases as the square of the radius, equation 
(1 - 49) may be transformed into a factor 

n B = (1-^!^* sin^ ) 2 (1 - 50) 

For the ideal windmill ^ as was numerically calculated from 
(1 - 30 a) together with (1 - 30). 

The calculated power was the power as if the number of blades were 
infinite. This power may be corrected through multiplication with the 
factor ng. 

In propeller theory, the propeller efficiency is defined 

n = LL (1-51) 

where V is the speed of the airplane, T-is the thrust, ft is the angular 
speed of the propeller and Q is the necessary torque. 
The tip losses diminish both torque and thrust. Since the reduction of 
the torque is smaller than the reduction of the thrust, there is a 
small drop in the efficiency of the propeller. For propellers this 
efficiency loss is minimized by changing the form of the blade tips 
(3, 5); the tips are rounded off in an elliptic form. This results 
in even more reduction of torque and thrust, however in such a way 
that the efficiency is the highest possible. Now, in windmill per- 
formance calculations, the situation is quite different. 
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Defining the windmill efficiency as n = leads us on very dangerous 
roads: if a relative reduction of the torque corresponds to a larger 
relative reduction of the thrust we have to conclude that the efficiency 
is higher. Although this might be true, the results would be bad since 
we would actually have designed a windmill with less than maximum power 
for the sake of obtaining the highest efficiency. It is therefore that 
the proposed efficiency be not used. Instead a coefficient should be 
used that relates the output to the power that is in the wind. 

C p = 1~ = (1 - 20) 

ipV^A ipV^A 

From this difference it immediately becomes clear that the negative 
effect of a finite number of blades is felt much stronger in calculating 
windmill power coefficients than in calculating propeller efficiencies. 
Although reduction of torque is accompanied by a reduction of thrust, this 
thrust is of no interest at all. Reducing the chords near the tips will 
lead to even more losses of torque and power. 
Appendix III shows a schematic representation of the total energy 
conversion on an element at radius r in the form of a diagram as used 
in hydraulics. 

. Theoretical design for highest power coefficient . 

Given some airfoil, and with the airfoil its C Q - C L characteristic 
(5) , we may use fig. 1.7. 

As stated on page 11 and 15 the curves in this graph are best interpreted 
as collections of maxima of Cp-A characteristics ; for a certain Cg/C^-ratio 
we may choose from fig. 1.7 for which tipspeed ratio we will design a 
windmill. Neglecting for the moment the occurring tip losses, for maximum 
power we will choose the tipspeed ratio that shows the highest possible 
Cp value. 

From paragraph 1.1.2. we know that a relation exists between a and a' 
(1 - 22) and that, for an- optimum design, relation (1 - 23) must be 
realized between a and a 1 . 



(l-a)a = a' (l+a l ) 

a'= (l-3a)/(4a-l) 



{1 - 22) 
(1 - 23) 
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From fig. 1.6. a trigonometric relation was derived 



Having chosen the tipspeed ratio X , the value x p is known for every section 

of the blades (A r = xr/R). 

With relations (1 - 22), (1 - 23) and (1 - 30) it is possible to express 

both interference factors a and a' as well as x as a function of $ 

a * = - 05 * (1 - 52) 

cos <H1 



a ,* = 1-cos <f> (1 . 53) 

2cos 0-1 



x = sin <P (2cos $-1) (! _ 54 j 

r (l-cos*)(2cos(N-l) 



where, in order to indicate that expressions (1 - 52) and (1 - 53) are 
the desired values of a and a' for an ideal windmill, designed with 
theories based on the assumption of an infinite number of blades, the 
notations a* and a'* are introduced. 

Since X p is known these relations give the optimum values of a and a'. 
By equating thrust as derived from momentum theory to the thrust as 
derived from blade element theory, relation (1 - 35) was found. 

aC. COSd> 

* = _Jz I (l - 35) 

1_a 4. sin 2 * 

Combining (1 - 52) with (1 - 35) leads to 



C L = 4 ( l-cos4>) (1 - 55) 

Thus for a blade section at a known radius r and speed ratio x r the product 
of chord, number of blades and lift coefficient may be calculated. From 
the airfoil characteristics C L (a) and C L (C Q ) the value a is chosen and 
with equation (1 - 29) the blade angle 6 is found. 
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■ 



6 « * " a (1-29) 

A constant angle of attack a (with minimum C^/C^-ratio) results in an 
irregular form of the blade. By varying a and therewith C L> a small 
increase of C D /C L may be tolerated for realizing a blade form with 
constant or linearly varying chord. In all these design methods 
no consideration is given to tip losses and the expected output should 
therefore be corrected with the factor n B as discussed in the last 
paragraph. A second method to find this same power coefficient is by 
making use of equation (1 - 47) and (1 - 48) for calculating the flow 
conditions on the blade elements in the way described on page 13. 
From these conditions the torque contributions from the elements may be 
calculated with equation (1 - 11) expressed in the factors a and a'. 

$ = p (2nr)(l-aF)V oo 2a'F fir 2 (1 - 56) 

AQ 

aCq _ jpV^nRS ^ dC Q _ x dQ 
Ax Ar ~ dX 



R 



|pV 2 nR 2 dr 



Thus (1 - 56) may be written dimensionless 

^2 = 8(l-aF)a'FA x 3 (1 " 57) 

dx 

Note that a and a 1 are reduced by the factor F since these a and a' 

are not the average interference factors. 

See for F equation (1 - 40), (1 - 43) and (1 - 30a). 

By numerical or graphical integration Cg may be found and multi pi i eating 

Cq with A yields C p for one value of A . 

This process of iterations and integration must be performed for various 
a 's and the result is the C p -A characteristic of the windmill. 
The above described method for design (1 - 55) neglects the effect 
of the tip losses. Thinking about what the tip loss models say about the 
flow we came to the conclusion that to some degree the consequences of 
the specific phenomena, that occur at the tips, might be compensated 
for by an appropriate rotor desian. As earlier stated on oaae 18 the tip 
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speed correction physically means that the maximum change of 

axial velocity only occurs on the vortex sheets and that the average 

change of velocity is only a fraction F of this maximum change. 

If a change of axial velocity of 2aV oo is thought to occur, the average 

change will only be 2aFV oo . From the simplest momentum models we know that 

the maximum power extraction is possible when the axial velocity is 

reduced to one third of its original value. If we want to realize this 

reduction we should design in such a way that the average kinetic energy 

in the wake flow equals that of a flow with velocity equal to one third 

of the wind speed. Consequently the reduction on the vortex sheets must 

be higher than the desired average. 

Design formula (1 - 55) was derived from (1 - 35). This last relation 
was the result of equating thrust as found from momentum theory to thrust 
from blade element theory. If we want a larger reduction of the velocity 
on the vortex sheets, we must therefore enlarge the numerical value of 
factor a* (as calculated from (1 - 52)) in the blade element part of 
equation (1 - 35). Taking again equations (1-1) and (1 - 11) we may now 
write, denoting by a* the desired average interference factor as calculated 
in (1 - 52) 



V (l-a*)2a*V = V^q-aV 

ir 2 sin 2 ,}, C v 



(l-a*)a* 
(1-aV 



aC^COStJi 

4sin 2 <|> 



(1 - 58) 



Factor K must be applied to the velocity change on the vortex sheets to 
obtain the desired average velocity reduction of the flow. If average 
velocity reduction would coincide with average kinetic energy reduction, 
then we could write 1/F for K in (1 - 58). Since in the energy the velocity 
is involved as a square, a reduction of the velocity change leads to a 
relatively greater reduction of the energy. If, in (1 - 58), we write for 
K=l/F, the result will be that the average velocity in the wake is one 
third of the undisturbed flow velocity but the remaining kinetic energy 
is still too high. Therefore K should even be larger than 1/F. Although 
we know that K should be larger than 1/F, (1 - 58) is rewritten with 
K=l/F as a first approximation. 
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m (l-a*)a*r . 4sl to (1 - 59) 

L (F-a*) 2 



Factor F may be calculated with (1 - 41) and (1 - 43). For elements near 
the tip F approaches to zero. In (1 - 59) it may be seen that already at 
F - a 1/3) the desired value of aC L becomes infinite and if we really 
want the velocity change to be ZaV^ we must, for values of F<1/3 S feed 
power to the system. 

Cons us ions and remarks 

- To describe the tip losses it is possible to use different modifications 
of the Prandtl model . 

- The modification that seems the most realistic, is using the angle ^ 
of the flow in the slipstream to calculate the distance between the 
vortex sheets of the slip stream. 

- Reduction of efficiency losses due to tip losses is realized in propeller 
theory by diminishing the tip chords. In windmill design not the inverse 
propeller efficiency but the power coefficient should be maximized. This 
leads to relatively large chords near the tips of the blades and only for 
values of F>a total compensation is possible. 

- The analysis is based on two-dimensional airfoil theory. Large changes 
in chords along the radius lead to the situation that the model is not 

any longer adequate because the flow becomes essentially three- dimensional. 



-27- 



SELECTION AND DESIGN OF THE ROTORS 

Windmills may be built from several materials varying from very simple 
to very sophistocated. Examples are cotton, wood, bamboo, steel, dacron, 
epoxy etc. These materials may be used in various ways. The material, 
together with the way it is going to be used, demands a certain level 
of technology. Therefore, in (6), different materials and combinations 
of materials were suggested for experiments: 

1. thin arched steel plates (see fig. 2.1.) 

2. double steel plates 

3. masts + poles 

4. sail wings 

5. sail trousers 






a 



b 



c 






d 



e 



f 



fig. 2.1, Simple types of airfoils 



a. thin steel plate 

b. double steel plate 

c. mast + pole 



d. sail wing 

e. sail trouser 



f. modified sail trouser 



See (6) for more information on these airfoils and the there given 
references. 
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For the following reasons the information on windmills built from 
these airfoils is very limited: 

a. only incidentally measurements and experiments have been done 
and recorded 

b. information on lift and drag properties of these airfoils, 
from which behaviour and possibilities of the windmills might 
be predicted, is very limited and questionable. 

Therefore, as a start, comparative experiments were suggested. 

To do these experiments an open wind tunnel was designed and built. 

See chapter 3 for a description of thetunnel. 

The windmill to be tested is placed in a free air stream. The test 
rotor diameters are 1.5 m (see chapter 3). 

As selection critirion a high value of the power coefficient is chosen; 
as a secondair critirion a high speed. 

With the help of the C D /C L -ratios from the literature study (6) and 
fig. 1.7. the following design speeds were selected: 

mast + sail (C D /C L £ 0.2) x design- = 2 > 3 

steel plates and 

sail trousers (C D /C L % 0.1) x design = 2 5 3 s 5 

double steel and 

sail wings (C D /C L % 0.05) x design = 3 9 5 s 7 



The second parameter that was fixed for design was the lift coefficient. 
As stated in chapter 1 paragraph 3, a fixed as well as a variable C^-value 
may be chosen. In our case we fixed the value of on that point of the 
C L -cx-characteristic that was thought to correspond with minimum Cg/C L - 
ratio. ^6). This resulted in the following design values for C. : 
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mast + sail 




0.8 


steel plate 


C L = 


0.8 or 1.0 


double steel plate 




1.0 


sail trousers 


c L = 


1.0 


sail wing 




1.2 



With formulas (1 - 54) and (1 - 55) it may be seen that, having fixed 
tip speed ratio A and lift coefficient C L , the value of the solidity a 
at every radius may be calculated. Now one more choice must be made 
because the solidity includes two variables: number of blades B and 
chord c. For low design values of x the solidity a will be relatively 
large. Taking for these low A -values a low number of blades leads to 
very large chords of the wing and has a negative influence on the maximum 
power coefficient (see chapter 1). Taking a high number of blades for 
designs with high value of A leads to structural problems because the 
blades become too thin. With formulas (1 - 50) and (1 - 54) it may be 
seen that the influence of the number of blades on the power coefficient 
is relatively small for high values of A. For these reasons the number 
of blades is in general lower for higher A-values. Weighing the above 
mentioned effects, several rotors with different numbers of blades are 
proposed for experiments. A total picture of the proposed test program is 
given in table 2.1 on the following page. 
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mast + sail 



steel plate 



double steel plate 



sail trousers 



sail wing 



c. 




8 










A 




2 

Cm 




A 


» 


6 


A 






R = 
u 




> 




C L 




0.8 










X 




2 


B = 


4 






X 


= 


3 


B - 


3 


9 


4 


A 






R - 


? 


9 


3 


r 

L l 




i n 

1 « u 










A 






R = 

D 


c 






C L 




1.0 










\ 
A 






R - 
u 


3 






\ 
A 




5 




2 






X 




1 


B = 


2 






r 

\ 




1 










A 




3 


B = 

u 


2 


» 


3 


A 




v) 


B 2 








C L 




1.2 










X 




3 


B = 


2 




3 


X 




5 


B * 


2 






X 




7 


B = 


2 







Table 2.1. Proposed test program. 



-31- 



During the period in which the test rotors were to be selected and 
designed the effects of the number of blades was not understood in 
the way explained in paragraph 1.3. and in the design formula the 
number of blades was not taken into account. Therefore only basic 
formulas (1 - 54), (1 - 55) and (1 - 29) were used in designing the 
test rotors. 



x = sin<j)(2cos<fr - 1) (1 . 54 j 

r (l-cos<f))(2cos<f> + 1) 



o C L = 4 (l-coso); a = |^ (1-55) 



= <f> - a (1-29) 



The result of evaluating these formulas for one combination of radius 
R, design-x , design-C^ and number of blades B is that, for every radius 
r on the blade, the blade setting $ and the value of c is known. The 
value of c is the length of the chord in the two-dimensional flow model 
of the blade section. \1) suggests that this chord should be set out 
perpendicularly on radius r. The way the air flows on the blades depends 
on the speed ratio of the element under consideration. For values of x 
near zero the path of the air is largely determined by the wind speed, 
while for high values of x the path becomes totally determined by the 
circumferential speed of the blade element. For low values of x the path 
thus becomes a straight line while for high values of x the path will be 
circular. Although the flow is circular we do not have to give up the 
two-dimensional model: just imagine that the plane of the intersection 
in two-dimensional modeling is now curved with a curvature depending 
on the speed ratio of the element. 

Fig. 2.2. shows the difference between setting out the chord as a 
straight or as a circular line. Point 0,0 is the center of the rotor. 
The point that is found with the straight line is called xl , yl and 
the point from the circular line is called x2, y2. A compromise weighing 
the speed ratio of the element at radius r may be found with the following 
formulas, derived from linear interpolation. 



x3 = x2 + (xl - x2)/(l + |x r ) 



(2 - 1) 



y3 = y2 + (yl - y2)/(l + f X p ) 



(2 - 2) 



y 4 



x3,v3 



_f----f xl.yl 
x2,y2 ^ 

V 
W 
\\ 

V 
Vi 
V 



Fig. 2.2. Setting out the blade chord. 



In Appendix IV results are shown for several designs. The effect of the 
"path correction" is highest for low design values of A . The effects may 
be neglected for elements withA r >2. 

The results of this analysis is not a significant increase in efficiency 
but a reduction of the blade chord at those stations where the chord is 
largest. This will result in less material use and in a more simple 
manufacturing. In this way the general plan form of the blades of the 
windmill design, according to table 2.1., were derived. Appendix V 
shows some examples of these designs. Note that they are calculated 
with constant C L ~ value and that the form may, if needed.be linearized 
as stated in paragraph 1.3. 

Of the five basic airfoil types, as mentioned on page 26, until now 
only two have been tested. The six steel plate windmills as suggested 
in table 2.1. were constructed and tested. A start was made with the 
sail trousers; two two-bladed and one four-bladed wi ndmi 1 Is were constructed 
and tested. 



2.1. Construction of the steel plate rotors 



From r ,4) and {7) it was concluded that the optimum arching of a steel 

plate, with respect to C D /C L -ratio, is about 7%. With the chord c, 

as calculated at r = 0.7 R, the radius R. that corresponds with the 

a 

calculated c and the desired value of the arching (7%), was calculated. 
See fig. 2.3. and table 2.2. 



c 




*- R = 1.82 x c (r = 0.7R) 
a 



Fig. 2.3. Calculation of the arching radius. 



Steel plate (thickness = 1.5 mm) is rolled until a cilindrical form 

with radius R a is res ized. With A ^ = £ xx and equation (1 - 54) the 
a r ft 

angle <j> at radius r = 0.2 R and at r = R are calculated. For a design 
with a constant C^value, the angle of attack a is for every radius the 
same. This means that the twist, which is the difference in blade angle 
at different radii, must be taken equal to the difference in angle <f>. 
Note that the smallest r (=0.2R) is arbitrarily taken as a measure for 
the hub of the rotor. The twist in fact is non-linear but since it is 
difficult to produce such a nonlinear twist an approximation is realized 
in the following way. A paper model of the blade form is put on the 
cilinder in such a way that the axis of the cilinder makes an angle with 
the axis of the paper blade. From the radius R a of the cilinder and 
the desired twist it may be calculated how the paper must be put on the 
cilinder. See fig. 2.4. When the form of the blade has been drawn in this 
way on the cilinder the paper may be removed and the blade may be cut 
out with a band-saw. In this way a blade has been produced with the 



right chords and a good approximation of the blade angles. Note that 
the radius of the cilinder was calculated from the chord at r = 0.7R; 
this means that for r<0.7R the arching will be stronger. See table 2.2. 
where the arching at r = 0.2R is given. To be able to change the blade- 
setting of the blade for testing, the blades were not directly welded 
on a hub, but were welded on footplates that were bolted on the hub. 
In this way only two hubs were needed: one for the two- .and four-bladed 
windmills and one for the three-bladed windmi lis. See fcr example the 
construction on the photo of Appendix VI. 




blade tip 



blade axis 



cilinder 
axis 



arched steel plate 



paper model of the blade form 



center of rotor 
The distance q may be calculated from 



q = 



R 

R-0.2R 



{■ 



(r=0.2R) v (r=R) 



71 R 

m ' a 



Fig. 2.4. Drawing of the blade form on the arched plate. 
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1.0 
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1.333 


7.54 
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17.14 


0.11 
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table 2.2. Design figures for steel plates 

2.2. Construction of the sail trouser rotors 

Was the information on arched steel plates limited, on the sail 
trousers there is no information at all. Therefore we decided to design 
the sail trouser rotors in such a way that many parameters (twist, blade 
angle, number of blades etc.) could, if neccessary, be varied. Since for 
two- and four bladed designs one hub could be used, the construction 
of a three- bladed rotor was postponed. On fig. 2.5. a sketch is shown to 
illustrate the way a sail trouser rotor was constructed. On the hub angle 
steels r £ or 4) were bolted; an angle steel forms the foot of a blade. 
On this angle steel a socket was welded. 

In this socket a tube can be fixed in such a way that later this tube may 
be rotated around its axis, for twist changes. At the other end of the 
tube, half of it was cut away over a length of approximately 10 cm. On 
this end of the tube a tip plate was mounted that has a width equal to 
the desired tip chord. The tip plate was bent at 25% of the chord 
from the leading edge, in such a way that, with the tube, a thickness- 
chord ratio resulted of 10%. From front and rears ide of the tip plate 
a leading edge and trailing edge were formed with steel cables connected 
with the foot plate. Over the construction of pipe and cables the trouser 
is drawn and fixed at the foot of the blade. The cable constructions 
only allow linearly varying chords and blade angles. Therefore the angles 
and chords, as calculated with constant C L over the blade, were treated 
in the following way. 



36 




g. 2.5. Sketch of the construction of a sailtrouser rotor blade mounted 
on a hub on the rotor axis. 
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At first the chords and blade angles were linearized between the chords 
and angles at radius r = 0.2R and r « R. The difference between the 
desired chord and the linearized chord is rather small, but the blade 
angles change so much that it may be expected that this poor approximation 
will result in bad performance. The outer parts of the swept area "catch" 
a relatively great part of the wind (the swept area is proportional to the 
square of the radius). The aerodynamic efficiency is low at small radii, 
because the speed ratio is low. A better result may be expected if the 
linearization of the blade angle is done between r = 0.5R and r = R. 
The change of the desired angle is small between 0.5R and R (see Appendix 
VII). This results in a small twist compared to the twist as calculated 
from linearization between 0.2R and R. The steel cable can be fixed at the 
foot plate at different places; in this way it is possible to change the 
average thickness of the blade. Two rotors were built in this way, both 
designed forx= 3 at = 1.0. One rotor has four and one has two blades. 
With two blades of the four bladed rotor it was possible to redesign 
the blade angles in such a way that another two bladed rotor was available 
for 4.5. The actual setting of the blades may be found in chapter 4 
together with the test results. 

2.3. Stress analysis 

The blades of a wind mill must be strong enough to resist two loads: 
centrifugal forces and bending moments. It will be shown in this paragraph 
that during operation the roots of the blades only have to resist the 
centrifugal forces and that the root tensions due to bending moments are 
relatively low. For when the rotor is not in operation, the bending moments 
must be calculated for storm conditions. 

2.3.1. Centrifugal forces 

The blade is modeled in the way indicated in fig. 2.6. The blade model has 
a constant cross-sectional shape. 

R 




Fig. 2.6. Model of centrifugal load. 
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The centrifugal force is then 



F c . m Z r =Pm A R .« =P m A4^ (2-3) 



-IT - "a —T- (2 " -») 



p = density of construction material 
m 

A = cross-sectional area 

ft = rotor angular speed 

o = stress at the root of the blade 

For a given material (wood, steel, etc.) the material properties p ffl and 
permissible are known andx max may be written as f(V«) 



2a. 



= 1W — EST!" (2 - 5) 

max 1/ p ( 



3 3 

For normal steel P|T] = 7.85 x 10 (kg/m ) 

and if °r, own is taken 160 * 10^ (Nm' 2 ) 
perm v ' 

the construction is safe against centrifugal loads if 

X <|00 (2-6) 



2. Bending moments during operation 

From chapter 1 we know the thrust on an element Ar 

AT - P (2nrAr) U (V^-Uy) (1 - 1) 

This indicates that the thrust is approximately a linear function of r. 
Then AT may be approximated by 
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The load on a blade under operating conditions may be modeled as 
shown in fig. 2.7a. Fig. 2.7b shows an element of this blade and the 
forces and moments acting upon it. 




fig. 2.7a. Load on blade under operating conditions 



AT T 2f 



F 




fig. 2.7b. Forces and moments acting upon a blade element Ar at radius 
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AM _ T F Ay _ AT aFc Ay 

AT " " 7 + C 2 AT T ~T" AT 



dF = - T 2 + F c 2 y 



t /f\ r a Fm /n 2 2 % 



With Pbl =^1»* 



F becomes 
c 2 



R 
F r = / Pbl fi 2 rAr = (R 2 -^) 
c 2 r * 



Combining (2-10), (2-11) and (2-13) gives 

dM - - Fm (R 2 r 2 \ + Pb1 ^ fR 2 r 2 \ u- 



The differential equation for elastic bending 



d y „ M 

2 ~ y ft 

dr^ Li 



differentiating (2-15) gives 



dM _ F r ... 

dF - E 1 y 
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Equating from (2-14) to (2-16) results in a third order differential 
equation with not constant coefficients 

EIY" ' - B^L (R 2 - r 2) y' = - 1$ (R 2 -r 2 ) (2-17) 



writing 



y = I (2-18) 



and X = £ (2-19) 

with ^ = ( )' (2-20) 

equation (2-17) becomes 

- Phl fi2R4 9 F m R3 o 

HrET- ^ " X > Y ' = " IFF (1 " X > ( 2 " 21 > 

? 4 
P hl O 

with Cj = D ^ EI (2-22) 



F R 3 „ 

and C x = (l-X^) (2-23) 



Y" 1 - Cj (1-X 2 ) Y' = - C 2 (1-X 2 ) (2-24) 



Conditions : 

atX=0:Y=0,Y'=0 (2-25) 
at X = 1 : M = + Y" = (2-26) 

For this differential equation no analytical solution is available and 
until now no numerical evaluation has been made by us. 
When from equations (1-24) - (1-26) Y" (0) is found the bending moment 
at the root may be calculated with (2-15) and the root stresses may be 
calculated with 
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M 



(1-27) 



Where W b is the resisting moment of the root section. 

In (9) it is stated that during operation the centrifugal forces will 
offset the bending loads on the blade and that the centrifugal force 
as dealt with in paragraph 2.3.1, is the main load. From figure 2.7 it 
can be seen that this offsetting depends on the angular speed. It is 
suggested that the above mentioned numerical evaluation be made to find 
out below which design tipspeed ratio the centrifugal forces. It is 
expected however, that no problems will arise. 



2.3.3. Bending moments for storm conditions 

During a storm, when the windmill is not operating, the blades must 
resist wind speeds of, for example, 40 m/s. The thrust on a blade 
element then is 



AT = 




D 



(2-28) 



where C Q is the drag coefficient. 
The bending moment becomes 



R 



M 



= J pV « 2 f C D crdr 



(2-29) 
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As an example the following figures are given 
\ =10, B - 1, c = 0.1 m, r _ 0.75 m c and f> n may be taken constant 
over the blade span, and the value of C n is 2. Then M becomes 



M 



= I x 1.3 x 1600 x 2 x 0.1 x \ 0.75 2 = 58.5 Nm 



With equation {2 - 27) the bending stresses, for storm conditions, 

may be calculated 

operating: 

M M 6 x 58.5 , c - 1rt 7 M -2 1C - A „ -2 



° b = \ " JTh 2 " ~ 0.1xl.5 2 xl0" 6 



= 156 x 10' Nm u = 1560 Nmm 



P 9 9 9 

not operating: T c = ^ x v f s 628 Nmm " 

This enormous root stress shows the importance of the protections of 
windmills against storms. The advantage of constructions built with 
sails is obvious: the sails may be taken off or reefed in such a way 
that the blade area is reduced. The calculations above were made for 
a hypothetic device that was designed for a high speed as x =10. 
To get an impression of the stresses that must be expected for relative- 
ly slow running windmills, another example is given. One of the rotors 
that we built and tested has the following properties: 

B = 4 

X = 2 

c = 340 mm (=average chord) 

c = 280 mm (=chord at root) 

h = 1.5 mm (thickness) 

P m = 7850 kgnf 3 

V = 15 ms"* (operating conditions) 
40 ms * (storm conditions) 



R = 0.75 m 

-1 



E = 0.21 * 10 12 Nm" 2 



p d1 $ 4 kgm 
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Centrifugal forces 

^2,2 ( 2 " 4 > 



at A = 2 o c = 3.5 Nmnf 2 

-2 

at A =3 o c =8 Nmm 
Centrifugal forces for operating in a storm 



a c - X 2 V„ 2 (2 - 4) 

at X = 2 a =25,1 Nmm 

C 

_? 

at A = 3 a =56.5 Nmm 
u c 



Bending forces for storm conditions when the windmill is not operating 



R 

M = \ p V m 2 / C Q crdr (2 - 29) 





M * 7 ^-Vaverase * 7 r2 



= | x 1.3 x 1600 x 2 x 0.34 x | x 0.75 2 = 199 Nm 



a - M = M = 6 x 199 s li9 x i 9 Nm" 2 = 1900 Nmm" 2 

b ¥ b |ch 2 0.28xl.5 Z xl0" b 
o 



3.4. Cable tension calculations 



Like in the previous analysis, the load on the blade of a sail trouser 
rotor is supposed to vary 1 i nearly with the radius. See fig. 2. 8. a. and 
equation (2 - 7) 
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From chapter 1 we know that the maximum energy to be extracted is 
P = Cp^pV^im 2 = I^pV.W (2-27) 



This is possible if the axial speed of the flow through the rotor was 



U = (l-a)V m = | V_ (2 - 28) 



The thrust to be expected is therefore somethingin the order of 

1 V « 



and for an element Ar per blade 
2 

AT = 4^ 2nrAr ( 2 - 30) 

If the load is supposed to be carried fully by the leading edge cable, 
the deflection and tension S of the cable may be calculated in the 
following way. See fig. 2.8. a. and b. 




fig. 2.8. Cable tension model. 
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V 2 

AT = Sa$ = £ ~- g- 2HrAr (2 - 34) 



d* . 8npV ^ r (2 - 35) 

~3r ~ 9BS r 

* - !^ r 2 + C = C x r 2 + C < 2 " 36 > 



For small values of y (see fig. 2. 8. a) tan<fr % sin<}> % 
(2 - 36) may now be written 

{j£ = cV t Cl (2 - .37) 

integration of (2 - 37) gives 



y = | C*r 3 + ^r + C 2 (2 - 33) 



Conditions: r = y = 
r = R y = 

Combining (2 - 38) with (2 - 39) yields 



(2 - 39) 



C x , 3 n 2 



y = j-(r -R"r) (2 - 40) 



The maximum deflection occurs if y*= 0. This happens when r becomes 
R/vT. The maximum deflection is then 



8n P V 2 R 3 

^max = — ~ — (2 - 41) 

81V3 BS 
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If for y__ v we take arbitrarily a value of 1 percent of the blade span, 
the needed tension 5 for the leading edge cable may be caculated with 



S = 



100 x 




(2 - 42) 



81V? B 



In our sail trouser constructions, the tube, as described in 2.2., is 
placed at 25 percent of the chord from the leading edge. If a balance 
of bending moments is desired, the trailing edge cabletension should be 
taken one third of the value of the leading edge cable tension. 
With the desired cable tension the diameters of the cables may be chosen. 
In this way the cable diameter for the 4-bladed windmill was chosen 2.5 mm 
and for the 2-bladed windmill 3.5 mm. 

2.3.5. Tube torsion 

Suppose that the total force on the trouser blade construction acts on 
the tip of the blade at the leading edge and the that the thrust, as 
may be calculated by integrating (2 - 33), is representative for the load 
on the blade. The torque acting on the tube may then be calculated with 




. 2nrdr 




(2 - 43) 







With the following set of equations the torsion ip and the shearing 
stress t may be calculated 



M.R 



"t" 
GI 

P 



(2 - 44) 



T 




(2 - 45) 
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where 



G = 0.385 E (2 - 46) 

V 3!T < d outSn 4 > < 2 " 47 > 



V if ( d out 3 " d in 3 ) ( 2 " 48 > 



To illustrate these calculations results are given for the two bladed 
sail trouser rotor. 



Given: c. . = 0.06 m 
tip 

-3 

p =1.3 kgm 
V = 15 ms" 1 

CO 

R = 0.75 m 

B = 2 

d. = 21 x 10" 3 m 
in 

d out = 21 - 8 * 10 " 3m 

E = 0.21 x 10 12 Nnf 2 



the torsion $ and the sheering stress x become 



= 0.3° 



t = 8 Nmm" 2 
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Conclusions and remarks 

- A more accurate interpretation of the calculated chords is qiven. 

This does not result in a significant increase of the power coefficient 
but leads to less material use and a more s impel manufacture 

- A simple method for manufacture of twisted, arched steel plates is 
given. Six rotors were built of blades that were manufactured with 
this method. 

- First designs of sail trouser blades are made and two rotors were 
made according to these designs. 

- Calculation of stresses in windmill rotors is complex even if 

a static analysis were sufficient. For static loads no strength 
problems are expected. It is thought however that, in connection 
with expected fatigue problems, also an analysis with dynamic loads 
should be made. This analysis was not possible within the scope of 
this work. 

- The static load on the blades under normal conditions is low but a 
storm will destroy each windmill that is not protected in some way. 

- Stress calculations for the components of the sail trouser frame are 
simple and no strength problems are expected. 
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3. WIND TUNNEL AND MEASURING INSTRUMENTATION 

In this chapter the test unit and the objective of several parts of 
the used instrumentation are explained. The used instrumentation is 
specified in Appendix VIII. At the end of this chapter a survey is 
given of the necessary relations for calculating the desired quanti- 
ties A, Cp and Cq from the measured quantities. 

3.1. The wind tunnel 

To test the rotors, a wind tunnel was built as is schematically shown 
in fig. 3.1. A 2.2 m diameter, 45 kW fan is driven bij a direct 
current motor. The power for this motor is generated by a socalled 




/?////// 



Fig. 3.1. Wind tunnel 
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Ward-Leonard set. In this way a stepless adjustable fan speed is 
obtained and the flow velocity may be stepless adjusted between 
zero and maximum. As indicated in fig. 3.1. the air passes several 
objects before leaving the tunnel. The objective of these vanes, 
gauzes and pipes is to eliminate rotational movements and velocity 
differences in the air. Thus a free airstream is generated of 2.2 m 
diameter. The velocity range is from to 15 meters per second. 
The velocity variations over the cross sectional area, in the free 
airstream after the tunnel, are less than two percent of the average 
velocity and the flow is entirely axial. In this airstream the wind- 
mill is placed at 2 m distance from the tunnel opening. Tests have shown 
that at this distance from the opening only a few centimeters of the 
boundery between moving and not -moving air have interfered with each 
other and that the velocity distribution of the flow that passes the 
rotor area is still the same as in the tunnel opening. Other measurements 
have shown that, when the windmill is operating, the velocity distribution 
in the tunnel -opening is not affected by the presence of the windmill 
2 m down stream. It is therefore assumed that the flow is sufficiently 
similar to the wind in open air. 

To prove that the finite dimensions of the free airstream have 
no influence on the test results, the following experiment was done. 
A simple four bladed steel plate rotor was built with diameter equal 
to half the diameter of the windmills to be tested. This rotor was tested 
at 2 m distance from the opening of the tunnel. Then a tube, with diameter 
equal to half the wind tunnel diameter, was placed in the wind tunnel 
opening in such a way that an airstream was generated with similar flow 
conditions but half the diameter of the original airstream. In this 
small diameter flow once more the 0.75 m diameter test rotor was tested 
under the same flow conditions, as in the large diameter airstream, at 
a distance of 1 m from the opening. See fig. 3.2. 

As shown in fig. 3.3. the differences between the test results were not 
significant and it was concluded that the behaviour of the small rotor 
was not influenced by doubling the flow diameter. We therefore concluded 
that for the 1.5 m diameter rotors the same test results will be found 
in a 2.2 m as in a 4.4 m diameter free air stream. 
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The ratio between the rotor swept area of the 0.75 m diameter rotor 
and the cross sectional area of the free airstream of 2.2 m diameter 
is 0.12. This ratio is such that this situation will yield the same 
results as if the rotor were placed in a flow of infinite diameter. 
This means that, with respect to the rotor diameter / air stream 
diameter ratio, the test set up is representative for the atmospheric 
wind. 

A flow property that might affect the energy transformation by the 
rotor is the degree of tubulence of the flow. From the construction 
of the tunnel it is obvious that the flow is strictly not laminar. 
The occurring turbulence is of a rather small scale due to the fine 
wire gauzes. The natural winds are also rather turbulent but the 
degree of turbulence as well as its scale is different in different 
situations. Therefore no special arrangements were made to influence 
nor to measure the degree or the scale of the turbulence. From the 
blade chords of the rotors and the rotor speeds the Re-numbers for 
the flow on the blade sections at 0.7 R were calculated as a function 
of V . See table 3.1. For the used airfoil sections the boundary 

00 
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layer flow is estimated to be turbulent above Re = 10 . 
It may be expected that the turbulence of the free air stream from 
the tunnel has only a very small effect on the performance of the 
rotors, as V in the experiments was about 13 m s 



airfoil 


design 


B 


Re(r=0.7R) 


steel plate 


2 


4 


3.5xl0 4 xV 

00 


ii 


3 


4 


2.6xl0 4 xV 

CO 




3 


3 


3.4xl0 4 xV 

00 


ii 


5 


3 


2.2xl0 4 xV 

00 


ii 


5 


2 


3.2xl0 4 xV 

oo 


sail trouser 


3 


4 


2.0xl0 4 xV 

00 


ii 


3 


2 


4.1xl0 4 xV 

00 



Table 3.1. Re-numbers for several rotors as f(VJ at r=0.7R. 

Regulation and measurements of power output and input . 

As shown in fig. 3.1. the rotor is mounted on a shaft that is supported 
by a tube. This axis is connected with a D.C, generator by means of an 
elastic coupling. The field voltage of this generator may be adjusted 
stepless between and 240 = V while the resistive load that is 
applied may be adjusted between 0.27 n and 3.6 fl» In this way it is 
possible to choose any number of revolutions between zero and maximum 
for measurements and it is possible to find the total Cp-X characte- 
ristic. 

The generator is mounted in bearings in such a way that the generator 
as a whole is free to rotate around its axis except for the fact that 
a support arm prevents a rotation. The torque that is generated by the 
windmill is found by measuring the force that this arm of the generator 



excerts on a force transducer as shown in fig. 3.4. Note that by calcu- 
lating from this force the torque, the losses in the generator bearings 
and brushes are included. 



arm 




generator 



force — - 
transducer 




Fig. 3.4. Measuring the torque 

The force transducer is of the strain gauge type. The low voltage 
output of this transducer is amplified and then converted into a 
frequency. This is done with a voltage to frequency converter that 
generates 1000 pulses per second per volt. By counting these pulses 
during 10 seconds an integration of the signal is realized. In this 
way the average of the varying output signal is measured; the variations 
are caused by the generator acting together with the force transducer as 
a mass spring system. During the same 10 seconds that the torque signal 
is integrated, the number of revolutions is counted (6 pulses per 
revolution during 10 seconds gives the number of revolutions per minute). 
As explained in paragraph 3.1. the flow velocity may be adjusted by 
means of the field voltage of the D.C. current generator that feeds the 
motor of the fan. The power input for the windmill is 




(3 - 1) 



-56- 



Before the rotor tests were started the velocity field in the tunnel 
opening was measured and it was concluded that the highest local devia 



The velocity is measured with help of a Pi tot-tube; the pressure diffe 
rence from this Pi tot-tube is measured with a Betz -manometer. The 
Pi tot-tube is placed in the tunnel 2 m before the opening at a spot 
where the velocity of the flow is equal to the average velocity in the 
tunnel opening. To be able to calculate the exact value of the air 
density the temperature and pressure of the air in the laboratory are 
measured and with these the air density is calculated. 

3.3. Formulas for calculation of the tipspeed ratio, power coefficient and 
torque coefficient from the measuring data . 

The pressure difference ap from the Pi tot-tube is read from the mano- 
meter in mm water column. If h is the height of this column then 



tion was less than two percent of the average velocity. 




(3 - 2) 



p 



air 



p air ( t=0 ° C ' p=760 mm Hg) * * 76ff 



(3 - 3) 



where T 



absolute temperature and p = air pressure in mm Hg 



The tipspeed ratio X is 



nR 2nnR 



T~ = 60V 

CO I 



and the torque Q is 



Q = F x 1 = <5 x ANP x 1 



(3 - 5) 
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where ANP = number of counted pulses 

1 = length of arm on generator that excerts a force on the 
force transducer. 

s = coefficient found from calibrating the force transducer 
including amplifier and voltage to frequency converter. 

F = force as measured by the force transducer. See fig. 3.4. 
y 

Cp and Cq may be calculated with 



2nnQ 
T 11 



Calculation of the error in the power coefficient 



C P = 575p7 = lp(M)3/Z nR 2 (3 - 8} 



Neglecting the error in p as calculated by (3 - 3) and the error in R, 
we may write 



Fx! xn 

C p = Cx(-3 ) (3 - 9) 

(.P) 3/2 



The error in factor 5 in equation (3 - 5) is very small and will be 
neglected, 

C p = Cx(aNP x 1 x n x Ap" 3/2 ) (3 - 10) 
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dC, 




dl + 



dn + Sap dAp 



(3 - 11) 



Differentiating C p in the way indicated in (3 - 11) and deviding 
equation (3 - 11) by C p as given in (3 - 10) gives 



ou P a(aN 
7^T = ~W 



+ 



Al 
1" 



An . 3 Aj[Apl 
n 7 Ap 



(3 - 12) 



ANP is the difference between the pulses as counted at some speed where 
we want to measure the torque and the pulses as counted when the rotor 
torque is zero. The error in aNP is estimated in the following way. 
The arm of the generator was lifted 10 times from and placed on the 
force transducer; because the generator is out of balance, a force 
acts always on the transducer, even if no torque is applied to the 
generator. 

In this way 10 times the force on the transducer was measured by coun- 
ting during 10 seconds the pulses as generated by the voltage to 
frequency converter. Of the results of these 10 measurements, the standard 
deviation s was calculated. As maximum error in NP 2s was used. Before 
every experiment the number of pulses NP must be counted and the number 
of pulses as counted for one torque measuring must be diminished with 
this initial number. Therefore the error in aNP might be 4s. The experi- 
ment as explained above resulted in s=40. Thus a(aNP) = + 160. 
Al is estimated 1 mm 

An is maximum 1 rev/mi n 

a(ap) is estimated 0.2 mm water column 

One of the rotors was tested and the following figures were measured. 

Ap = 10.7 mm water column 
ANP= 49000 pulses 
n = 444 rev/mi n 
1 = 300 mm 
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Uith (3 - 12) the relative error in Cp can now be calculated 

= T50W + m + m + ! TCT7 ■ 0.0032+0.0033+0.0025+0.028-0.037 

Cp as calculated from the above given figures with the equations in 
paragraph 3.3. is 0.36. With the relative error of 0.037 this result is 
to be interpreted as 

C p = 0.36 + 0.013 

From the above given analysis it is concluded that the error in Cp is 
mainly due to the error in measuring Ap of the Pitot-tube. The influence 
of this error is about ten times higher than the error in each of the 
other ones and causes 75 percent of the possible error in Cp. Since all 
our measurements until now are done at the same flow velocity, the 
maximum Cp values of the other rotors will also include a maximum 
relative error of about 0.037. 

Conclusions and remarks 

- A test unit has been built for testing rotor models with diameters 
of 1.5 m. 

- Experiments have shown that, with the chosen rotor-airstream diameter 
ratio, the finite dimensions of the free airstream have no significant 
influence on the test results. 

- It is expected that the test results are representative for the 
behaviour of the same rotors in an atmospheric wind except for two 
aspects: degree of turbulence and the constancy of the flow. Compara- 
tive experiments should be carried out as soon as possible. 

- Re-number of the flow on the blades, calculated with the chord of the 

-1 5 
blades at r=0.7R at V ra =13 ms , vary for optimum x between 2.6x10 and 

5.3xl0 5 . 

- Flow velocity in the test unit may be adjusted between zero and 15 ms" 1 
and the load on the windmill may be adjusted in such a way that every 
point of the Cp-x curve of the rotor can be tested. 

- Maximum relative error in the maximum C D values will be about 3.7 percent. 
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4. ROTOR TESTS 

Various rotors were tested with several configurations and blade- 
settings. The objective of every test was to find experimentally the 
Cp-A and Cq-a characteristics of the rotor. These characteristics 
present the properties of the rotors in a form that allows comparison 
of these properties. 

4.1. Steel plate rotors 

In chapter 2 the selection and design of these rotors was motivated 
and explained. In (6) a C^-a characteristic of an arched steel plate 
is shown. Though this arched steel plate is not symmetric, it is. thought 
that with this information a first choice of the angle of incidence a 
can be made. For a non- symmetrical arched plate the value of a is 3° 
for C L = 0.8 and 4° for C,_ = 1.0. 

If the angle <t> between the relative velocity and the plane of the 
rotor were known then it would be possible with 

B = <J> - a (1-29) 

to calculate the desired blade setting 3. The angle <f> is only known 
for the ideal flow and may, for optimum energy extraction, be calcu- 
lated with (1-54). The result is shown in Appendix VII where <fc ideal 
is shown as f(A r ). Using these values for <j> a design value for b is 
found for the blade element at r = 0.7R. Table 4.1. shows the chosen 
values of s. 

The first rotor test showed a maximum Cp value at a smaller blade 
setting than the design value. We therefore concluded that the design 
value of a should be taken a little higher. For rotor 2 and 3 the angle 
of incidence a was taken 4°. This proved to be incorrect as will be seen 
from the test results of rotor 2 and 3. 
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Table 4.1. Determination of the blade setting 

The rotors were tested in the following way. The blades of a rotor 
were set at the v< >ut* of 6q from table 4.1. For about ten values of the 
number of revolutions between zero and maximum, measurements were done 
and Cp S Cg and x were calculated. Then the blade setting 6 was changed 
a few degrees in an arbitrary direction and again the rotor was tested. 
If the maximum value of C p proved to be higher, the change of 6 was 
repeated in the same direction, else e was changed in the opposite 
direction. In this way a rotor was tested for several blade settings, 
including the optimum one. This was done with each of the six rotors. 
The results are shown on the next pages in fig. 4.1. up to and inclu- 
ding 4.6. 

Of the rotors that are designed for x = 5 the leading edges were 
rounded off and the trailing edges were sharpened. Fig. 4.7 shows 
the effect of this measure on one of the rotors. 

The marked points in fig. 4.1 - 4.7 are actual measuring results. 
Different marks indicate different rotor tests as specified in the 
figures. Fig. 4.8 shows the collection of the C p -X curves of the six 
steel plate rotors at their optimum blade settings. 



3 (3(r=0.7R) = ll° 
-13° 



steel plate rotor 

\-design=3, B=4 ( X 16 

C -design =0.8 




Figure 4.2 



n Yes 

leading edge round , trailing edge sharp ? 

a No 

steel plate rotor 
X-design=5» B=5» 
C T -design = 0.8 
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Sail trouser rotors 

The design of the form of the sail trouser blades was described in 
paragraph 2.1. The following blades are available: 

- 4 blades, thickness 10 percent, maximum thickness at 25 percent 
of chord from leading edge; chords are linearized and are given 
by the formula: c = 0.25 - 0.19r 

Twist and blade setting are stepless adjustable. Average thickness 
can be influenced. 

These four blades will be referred to as ST-1. 

- 2 blades, thickness 10 percent, maximum thickness at 25 percent 
of chord from leading edge, chords are linearized and are given 
by the formula: = 0.456 - 0,338r. 

Twist and blade setting are stepless adjustable. Average thickness 
can be influenced. 

These two blades will be referred to as ST-2. 

With ST-1 a four bladed rotor (x=3) and a two bladed rotor (A =4.5) 
were built and tested. 

As explained in paragraph 2.1 the twist is linearized between 
r = 0.5R and r = R. At 5 places between 0.5R andR the angle cf> was 
calculated (see Appendix VII) and with linear regression the best 
straight line through these points (r,<jj) was determined. 
From this linear function $ = constant x r, the difference between 
the blade angle of an element at the blade foot and of an element 
at r=0.72 m was calculated; this difference is called the twist. It 
is estimated that the desired lift coefficient will be realized by 
the following angle of incidence 

C L = 1.0 at a = 4° 

The blade setting 3 may now be calculated with equation (1 - 29). 



3 = <f> - a 



(1 " 29) 
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In this way the desired blade settings and twists were calculated 
and the results are given in table 4.2. The rotors are numbered 
7, 8 and 9; numbers 1 until 6 were used for the steel plate rotors 
of paragraph 4.1. 





A 


B 


blade type 


twist(foot-0.72 m) 


e(foot) 


7 


3 


4 


ST-1 


15.0° 


24° 


8 


3 


2 


ST-2 


15.0° 


24° 


9 


4.5 


2 


ST-1 


11.4° 


16° 



Table 4.2. Design settings of the blades of the sail trouser rotors. 

The sail trouser rotors were tested in the following way. 
The four bladed design for x = 3 was set according to table 4.2. and 
tested in the same way as the steel plate rotors: at different 
points between zero and maximum number of revolutions, the necessary 
quantities were measured and C p , Cg and x were calculated. Then the 
blade setting was changed by 2 degrees and again the rotor was measured. 
In this way the curves in fig. 4.9. are obtained. Then the twist was 
changed by a few degrees while the blade angle at the tip was kept 
the same as the tip blade angle that corresponds with the best C p -x 
curve of fig. 4.9. 

Again the rotor was measured and results are shown in fig. 4.10. 
Twist changes only resulted in a small drop of C p and therefore the 
second design (ST-2-blades , x=3, B=2) was only tested for several 
blade settings with constant twist. Results are shown in fig. 4.11. 
The third design (St-l-blades , A =4.5, B=2) was tested the same way 
and results are shown in fig. 4.12. 



1 p(r=0.12R)=?2° 
2 4 ° 
26° 

Sail trouser rotor, type ST-1, 

\-design = B = 4 , -< 

C^-design =1.0 

twist(r=0, 12R-0.96R) = 15° 




3 4 

Figure 4-9 



tvist(r-0.12R-O.96R)= M7° 

■•15° 
-1?° 

Sail trouser rotor, type ST-1, 
\-design = 3, B = 4» 
C^-design =1.0, 
p(r=0.96R) = 9° 
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Interpretatien and discussion of the test results 

Table 4.3. shows a survey of the test results. Rotor numbers refer to 
the numbr^s as used in table 4.1. and 4.2. 



rotor number 


1 


2 


3 


4 


5 


6 


7 


8 


9 


x-design 


1 


3 


3 


5 


5 


5 


3 


3 


4.5 


6-design 


21 


13 


13 


7 


7 


7 


24 


24 


16 


Cp-maximum 


0.36 


0.37 


0.32 


0.32 


0.34 


0.37 


0.26 


0.22 


0.17 


\ — nn - *" "1 mum 

A ULs ~> I IIIUIM 


7 5 


*3 


o 


A 
*+ 


A 




7 A 


c. , X 


9 ft 


r-od ti mum 


18 


16 


13 


X X 


] 1 
X X 


Q 
7 




77 


1 7 


r at 
Pmax dt 


0.35 


0.33 


32 


30 


31 


32 


26 


20 


n u 
\j ■ x j 


8-design 




















x optimum 
at $ -design 


2.4 


3 


3 


4 


4 


4.1 


2.4 


2.1 


2.7 


C Q(x=0) at 
B-optimum 


0.160 


0.080 


0.057 


0.025 


0.021 


0.014 


0.043 


0.037 


0.009 


Ht' rt 


0.144 


0.123 


0.107 


0.080 


0.085 


0.082 


0.108 


0.105 


0.060 


6-optimum 





















1>jle 4.3. Survey of the test results of steel plate rotors 

number (1-6) and of sail trouser rotors number (7-9) 



Comparison of the values of x-design with the optimum values of X 
shows a remarkable effect: for the rotor designed for x=2 the optimum 
value of A is higher than the design value, for x-design = 3 the 
optimum value is equal to the design value, while for the other rotors 
the optimum value is lower. An explanation may be found in fig. 1.7. 
We interpreted the curves in fig. 1.7. as collections of maxima of C p -x 
curves that may be realized with fixed blade chords and blade settings. 
Although this interpretation is very helpful in selecting the design-*, 
a more accurate interpretation is that one ideal curve (for one specific 
C D /C L -ratio) is the collection of points of different rotors that have 
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m that point a tangent common with the tangent of the ideal curve. See 
fig. 4.13. 




Only if this tangent is zero the x-optimum will be equal to X-design. 
If x-design is not chosen atxwhere Cp-ideal (fig. 4.14) is maximum, 
then a tangent of zero is impossible for x = x-design, since this would 
imply that the design would have operating points outside the ideal 
Cp-x curve. This is not possible because the ideal curve gives the 
maximum possible Cp-values for one specific Cg/C^-ratio. 
In our design (x = 2, B = 4) we assumed C^/C^-OA. This assumption 
leads with fig. 1.7. to a choice of x-design = 2. The C Q /C L -ratio of 
0.1 is too pessimistic, as will be seen further on in this paragraph. 
It will there be argumented that the true Cp/C^-ratio is probably closer 
to 0.05. For this value ideal Cp-x curves are drawn in fig. 4.13 and 
4. 14. The maximum value of this curve is reached between X = 2.5 and 
> = 3. This explains the differences in deviations of x-optimum from 
X-design: \-design = 2 leads to a rotor with an optimum value of x 




2 4 6 8 10 



Fig. 4.14. Operating points outside the ideal C p -curve are impossible. 

higher than two because the tangent is possitive. See fig. 4.1. For 
A-design = 3 this tangent is ^ and A-design = A-optimum while for 
A-design > 3 the tangent is negative and therefore A-optimum <X-design. 
Comparison of 3-design with 3-optimum shows differences that may be 
influenced by both the angle of incidence a and the angle $ between the 
relative velocity and the rotor plane. The only steel plate rotor that 
shows a lower value for B-optimum than 3-design is number 1. 
This may be understood from the fact that for this rotor A-optimum 
is higher than A-design. Appendix VII shows that for higher A-value 
the angle <f> is lower and with equation (1 - 29) we can see that 6 will 
indeed be lower. The difference in X-optimum between rotors with number 
2 and 3 is insignificant as may be seen from fig. 4.3.: the difference 
between results of the rotor at 6 = 15° and e = 13° is smaller than the 
possible error as calculated in paragraph 3.4. 
Comparison of Cp-maximum and A-optimum at 3-design with Cp-max and 
A-optimum at B-optimum leads to the conclusion that the designed rotors 
were only slightly improved by correcting experimentally the blade 
setting. 
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For comparison of the measured Cp-values with theoretically possible 
Cp- values table 4.4. shows C p - values calculated for different C^/C^- 
ratios and with different models. The C p -values are calculated in the 
following way. For the optimum x of the rotor, Cp-ideal (Cq=0) is cal- 
culated with (1 - 24). This value is diminished with the result of 
(1 - 38) and the result is multiplied with the value rig from (1 - 50) 
Cp in table 4.4 indicates that i}g was calculated with the original 
model of Prandtl . See equation (1 - 50) and take sin<f>, = t~~T 4 

Clpml indicates that for this calculation the Prandtl model is modified 
by taking in equation (1 - 50) for <f>j the angle <j>. Cp m 2 indicates that 
(1 - 50) was used with angle ^ as explained under equation (1 - 50). 
When considering the figures in table 4.4. it is absolutely essential 
to keep in mind that 

- the measured values of Cp are measured at a real rotor, with a hub, 
with linearizations of for example the blade angles and with other 
concessions that have been made to the ideal design, 

- the maximum relative error in Cp is smaller than 4 percent 
(see paragraph 3.4. ) 

- the calculated Cp-values are values to be reached with a theoretical 
design that is only approximated by our real constructions. 

From the table we have drawn the following conclusions. 

- Only the second modification of the Prandtl model and the assumption 
that C D /C L = 0.05 result in a calculated C p -value that is higher than 
the measured Cp-value. 

- Results with the sail trousers suggest that the Cp/C^-ratio of the 
used sail trouser blades is higher than 0.1 but lower than 0.15. 

It is expected that the minimum C^/C^-ratio is 0.1, as was estimated 
for the design. 





steel 


pi ates 




C D /C L =0.05 


yC L «0.075 


c D /c L = 


.1 


design 


A-opt 


B C p 


-measured 








K Pr 


ml 


K m2 


c p C p 
r Pr ml 


C P 

m2 


2 


2.5 


4 


0.36 


0.348 


0.382 


0.419 


0.319 


0.351 


0.385 


0.291 0.320 


0.351 


3 


3 


4 


0.37 


0.361 


0.391 


0.422 


0.326 


0.352 


0.381 


0.290 0.314 


0.345 


3 


3 


3 


0.32 


0.332 


0.370 


0.411 


0.300 


0.334 


0.371 


0.267 0.298 0.331 


5 


4 


3 


0.32 


0.347 


0.377 


0.408 


0.300 


0.326 


0.353 


0.254 0.275 


0.298 


5 


4 


2 


0.34 


0.305 


0.347 


0.392 


0.264 


0.300 


0.339 


0.223 0.253 


0.286 


5 


4.5 


2 


0.37 


0.313 


0.348 0.388 


0.262 


0.294 


0.328 


0.214 0.240 


0.268 


Sail trousers 


c D /c L = 0.1 


c D /c L = 0. 


15 






3 


2.4 


4 


0.26 


0.291 


0.320 


0.351 


0.235 


258 


0.283 






3 


2.1 


2 


0.22 


0.188 0.246 


0.315 


0.169 


0.202 


0.238 






4.5 


2.8 


2 


0.17 


0.223 


0.266 


0.314 


0.141 


0.160 


0.181 







Table 4.4. Measured and calculated C n -values. 
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Conclusions and remarks 

- Comparison of test results with calculation results suggest that the 
Prandtl model is best used when modified in the way explained under 
equation (1 - 50) of chapter 1. 

- Cp/C^-ratio of the arched steel plates is probably 0.05 while Cp/C L - 
ratio of sail trousers is minimum 0.1. 

- Fig. 1.7. is a good help in selecting x-design values. The optimum X of 
a rotor designed for this selected value will only be equal to X-design 
if x-design corresponds with the value of X that shows, for a chosen 
Cp/C L -ratio, the highest C p -value in fig. 1.7. If a specific design 
speed is the first criterion for selecting a rotor, this speed is to 

be translated in x-design by the definition of X (X = fiR/V^). Though 
the optimum value of x will not necessarily be the same as X-design, 
selection of another value for X-design will not lead to a higher 
Cp-value at the design speed of the rotor. 

- With rotors built from simple materials high power coefficients are 
possible. 

- Deviations of blade angles from the ideal values result in small 
performance losses. 

- Only the four bladed design for X = 2 results in a rotor with a star- 
ting torque that is higher than the torque at maximum Cp. This property 
makes the design very suitable for constant torque loads: the same 
wind-speed that can start this rotor is enough for speeding up this 
rotor to the number of revolutions where Cp-maximum is reached. All 
other rotors would, under a constant torque load, need a higher wind 
speed for starting than is needed for working at X-optimum. 
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5. SUMMARY OF CONCLUSIONS AND RECOMMENDATIONS 

- To describe the tip losses it is possible to use different modifications 
of the Prandtl model . 

- The modification that seems the most realistic, is using the angle ^ 
of the flow in the slipstream to calculate the distance between the 
vortex sheets of the slipstream. 

- Reduction of efficiency losses due to tip losses is realized in pro- 
peller design by diminishing the tip chords. In windmill design not 
the inverse propeller efficiency but the power coefficient should be 
maximized. This leads to relatively large chords near the tips of the 
blades. 

- A more accurate interpretation of the calculated chords is given. This 
does not result in a significant increase of the power coefficient but 
leads to less material use and a more simple manufacture. 

- Calculation of stresses in windmill rotors is complex even if a 
static analysis were sufficient. For static loads no strength pro- 
blems are expected. It is thought however that, in connection with 
expected fatigue problems, also an analysis with dynamic loads should 
be made. This analysis was not possible within the scope of this work. 

- The static load on the blades under normal conditions is low but a 
storm will destroy each windmill that is not protected in some way. 

- Stress calculations for the components of the sailtrouser frame are 
simple and no strength problems are expected. 

- It is expected that the test results are representative for the beha- 
viour of the same rotors in an atmospheric wind except for two aspects: 
degree of turbulence and the constancy of the flow. Comparative experi- 
ments should be carried out as soon as possible. 

- Experiments have shown that, with the chosen rotor-ai rstream diameter 
ratio, the finite dimensions of the free airstream have no significant 
influence on the test result. 

- Re-number of the flow on the biases, calculated with the chord of the 

-1 5 
blades at r = 0.7R at V m = 13 ms , vary for optimum x between 2.6x10 

and 5.3xl0 5 . 

- Flow velocity in the test unit may be regulated between zero and 15 ms 
and the load on the windmill may be regulated in such a way that every 
point of the C p -x curve of the rotor can be tested. 
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- Maximum relative error in the maximum Cp value will be around 3.7 
percent. 

- Comparison of the test results with calculation results suggest that 
the Prandtl model is best used when modified in the way explained 
under equation (1-50) of chapter 1. 

- Cp/C^-ratio of the arched steel plates is probably 0.05 while Cp/C^- 
ratio of sail trousers is minimum 0.1. 

- Fig. 1.7. is s good help in selecting X-design values. The optimum X 
of a rotor designed for this selected value will only be equal to 
X-design if X-design corresponds with the value of X that shows, for 
a chosen C D /C L -ratio, the highest C p -value in fig. 1.7. If a specific 
design speed is the first criterion for selecting a rotor, this speed 

to be translated in X-design by the definition of X (X = ft R/V^) • Though 
the optimum value of X will not necessarily be the same as X-design, 
selection of another value for X-design will not lead to a higher 
Cp-value at the design speed of the rotor. 

- With rotors built from simple materials high power coefficients are 
possible. 

- Deviation of blade angles from the ideal values result in small per- 
formance losses. 

- Only the four bladed design for X = 2 results in a rotor with a start- 
ing torque that is higher than the torque at maximum Cp. This property 
makes the design very suitable for constant torque loads: the same wind- 
speed that can start this rotor is enough for speeding up that rotor 

to the number of revolutions where Cp-maximum is reached. All other 
rotors would, under a constant torque load, need a higher wind speed 
for starting than is needed for working at X-optimum. 
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! ff'-ct i-! iip Speed Ratio on the Induced \elociri<_ 
t'l-r f I mw \\ 1 th an i rrotat i ona I Kal.e . 
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APPENDIX III 



Schematic representation of energy conversion on blade element Ar at ra- 
dius r 



2tt 













V * \ 





ft 



d) = V . 2wrAr ; 

A = O.Spv/ ; 

_ P out _ n ■ Q = 



w 



0.5pV °A ; 



'a sm<J> L r 

, . Or . a = |c_ . 
r " IT ' ?Fr 5 

C D 

n w = (l - ^ cotcj>) ; 

tancf> = |p|ry x ; a = <|> - 3 ; C L = f (a) 

V* = (2TrrH2irrAr) , (2 „ rAr) . 

A s aer 

r 



C n = f (a) ; 
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APPENDIX IV 



i a 



i!el? oi stael pjLate blade 
B - 2, \-daalgn - 5, C L - 0.8 



straight path 
circular path 

compromise 



05 



15 



i 1 : 



.45 



.ss 



„ X t ... 

.65 .7 



! 

j- 

j Half of steel plate blade 

j- B - 3 , \-design - 5 , C..-0.8 



- 4b ■ ';> -aS 



-6 



Half of steel plate blade 
B - 4, \-deoign - 2» C.- 0.8 
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Half of symmetric steel pla^e bifida 
defligned for four bladed rotor, 
\-deaign - 2 , Cj-design - 0.0 



• 2S 



r>5 



Half of symmetric steel plate blade 
designed for a four bladed rotor, 
\-design - 3 . C -deals* - 0.8 




Half of symmetric steel plate blade 
designed for a four bladed rotor, 
\-design - 5 i C.-design - 1.0 
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APPBKDIX VII 



siaq?(2cosy-I ) 
(l-cos<p)(2cos9+l ) 



X 

r 


<¥> 


10.00 


3.807 


9*75 


3.904 


9.50 


4.006 


9.25 


4.113 


9.00 


4.227 


8,75 


4.347 


8.50 


4.473 


8.25 


4.607 


8.00 


4.750 


7.75 


4.902 


7.50 


5.063 


7.25 


5.256 


7.00 


5.420 


6.75 


5.618 


6.50 


5.831 


6.25 


6.060 


6.00 


6.308 


5.75 


6.577 


5.50 


6.870 


5.25 


7.190 


5.00 


7.540 


4.75 


7.926 


4.50 


8.353 


4.25 


8.827 


4.00 


9.357 


3.75 


9.954 


3.50 


10.630 


3.25 


11.402 


3.00 


12.290 


2.75 


13-322 


2.50 


14.534 


2.25 


15.975 


2.00 


17.710 


1.75 


19.830 


1.50 


22.460 


1.25 


25.773 


1.00 


30.000 


0.75 


35.420 


0.50 


42.290 


0.25 


50.642 


0.00 


60.000 



8 



r- 
1) 



Angle «p between relative velocity and plaae of the rotor versus 
the speed ratio of an element at radius r for a windmill with a 
flow equal to the flow of an ideal windmill. 
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List of used instruments* 

j 

- Electronic Counter ( 6 digits ), Hewlett/Packard, 5532A 

- Electronic Counter ( 4 digits ), Venner Electronics Ltd, TSA 6634 

- Static/Dynamic strain indicator, Pe.ekel N.V. , 540 DNH 

- Voltage to frequency converter, Dymec, DY-2211A/B 

- Force transducer, H.B.M., range 0-10 kg, Ul 

- Force transducer, Iweco, range 0-100 kg, ED335 

- Oscilloscope, Hewlett/Packard, 130C 

- Digital volt meter, Teletec Airtronic, TE 350 

- Betz micro manometer, van Essen 

- Pitot tube, Iweco 

- Photo-electric cell, v.d. Heem 

- Amplifier, v.d, Heem, 9818 

- Barometer, System Paulin A.B. 

- Laboratory thermometer. 



